Automated Chirp Detection with Diffusion Entropy: 
Apphcation to Infrasound from Sprites 
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We study the performance of three different methods to automatically detect a chirp in background 
noise, (f) The standard deviation detector uses the computation of the signal to noise ratio. (2) 
The spectral covariance detector is based on the recognition of the chirp in the spectrogram. (3) 
The CASSANDRA detector uses diffusion entropy analysis to detect periodic patterns in noise. All 
three detectors are applied to an infrasound recording for detecting chirps produced by sprites. The 
CASSANDRA detector provides the best trade off between the false alarm rate and the detection 
efficiency. 
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I. INTRODUCTION 

Chirps are periodic signals with an instantaneous fre- 
quency changing in time. Chirps are produced by a va- 
riety of sources: from lightning generated whistlers 
to the acoustic emission of bats ^ and whales 0. Re- 
cently the presence of chirps in infrasound recording has 
been associated with the occurrence of sprites over thun- 
derstorm clouds 3. Several methods of chirp detection 
have been developed, operating in both the time domain 
(multiple frequency tracker '5^ and recursive least square 
algorithm (Qj) and the frequency domain (Page's test 
and Hough transform 

In this work we introduce two new methods for chirp 
detection. The spectral covariance detector and the CAS- 
SANDRA detector. The spectral covariance detector 
operates in the frequency domain and uses the corre- 
lation between different frequency bins of the spectro- 
gram [l3 | as an indicator of the chirp occurrence. The 
CASSANDRA {Complex Analysis of Sequences via Scal- 
ing AND Randomness Assessment) detector operates in 
the time domain Ullil and uses diffusion entropy anal- 
ysis 0, We study the performance of these two 
detectors and compare the results with a standard signal 
to noise ratio detector (standard deviation detector). 

Sprites and other recently discovered Transient Lu- 
minous Events (TLEs) above thunderstorm clouds, such 
as elves blue jets an d g igantic jets 0| are 

the subject of intense research [l^. TLEs connect the 
lower layer of the atmosphere (troposphere: below 10 
km) where the weather activity occurs with the upper 
levels of the atmosphere (80-100 km). Knowledge of the 
sprite occurrence rate is of primary interest to address the 



'Correspondent Author. Email: eexmi@bath.ac.uk. 

Web Page: 

|http:/ / staff, bath, ac.uk/ eexmi / mywebpage /home / mywebpage.htm | 



relevance of these phenomena and their global impact on 
the atmosphere. 

The outline of this work is the following. In Section UTI 
we briefly discuss the infrasound signature of sprites. In 
Section imi we discuss the details of each detector and its 
performance in sprite detection. In Section llVI we draw 
our conclusions. 



II. SPRITES AND THEIR INFRASOUND 
SIGNATURE 

Sprites are TLEs with a typical duration from a few 
milliseconds up to a few hundred milliseconds. They are 
generated by the electric field pulse of a "parent" posi- 
tive cloud-to-ground (-I-CG) lightning discharge |23|. The 
vertical extension of sprites is ~45 km, starting from ~40 
km up to c:i85 km, while their horizontal extension can 
range from 20-50 km. Since the first optical observa- 
tions (15(1 , Sprites have been observed over thunderstorm 
clouds in North America [23, Europe and Japan 
|23'|. Electromagnetic signatures from sprites have been 
reported in the Extremely-Low Frequency (ELF) range 
(lOHz-SkHz) and with Earth-ionosphere cavity reso- 
nances :25|. 

The possibility that sprites could generate an infra- 
sound signature was first suggested by Lizska jl^ . The 
first report of sprite signature is by The sprite sig- 
natures are located in the 1-10 Hz frequency range and 
in many cases a linear chirp of increasing frequency with 
time is observed. This signature is caused by the spatial 
extent of the sprite (from 20 to 50 km) its orien- 
tation with respect to the infrasound station, and the 
reflectivity properties of the thermosphere 0| . Pressure 
waves generated from different regions of the sprite will 
be reflected at different altitudes in the thermosphere 
with different absorption and dispersion properties be- 
fore reaching the infrasound station. The net result is 
that pressure waves coming from the nearest end of the 
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sprite will arrive first at the station with a low frequency 
content. Pressure waves coming from the farthest end 
of the sprite will arrive later at the station with a high 
frequency content. Sprite signatures which show an im- 
pulsive feature instead of a chirp, are the result of a small 
spatial extension or of the alignment with the infrasound 
station (regardless of spatial extent). 

The data set used to test the automated chirp detec- 
tors is shown in Fig. ^ It is the original signal recorded 
with a 20 Hz sampling rate at the infrasound station in 
Flers (210 km West of Paris, France) from 2:30 UT to 
4:00 UT on the 21st of July 2003. During this hour and 
half an intense thunderstorm occurred in central France. 
Optical observations reported 28 sprites in the thunder- 
storm region. From these 28 sprites, 12 signatures were 
detected in the infrasound recording, 9 signatures were 
chirps with an average duration T~12-15 seconds as iden- 
tified by visual inspection. 

The sprite signatures are not visible in Fig. ^ because 
their intensity is very small compared to those of the 
slow pressure (/<1 Hz) fluctuations caused by the wind. 
But the spectrogram in the range 1 — 10 Hz as in Fig. [21 
shows the chirp signature of a sprite. Therefore before 
applying any detection method, we high pass filter the 
infrasound recording in order to eliminate the very slow 
wind fluctuations (/<1 Hz). 
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FIG. 1: The infrasound recording at the Flers station from 
2:30 to 4:00 UT on the 21st of July 2003. The values represent 
the pressure fluctuations (in Pascal) registered by a micro- 
barograph. 




FIG. 2: A chirp signature of a sprite is identified plotting 
the energy density (spectrogram) E{t, f) of the signal in the 
time-frequency plane for the 1-10 Hz frequency range. This 
plot refers to the 2 minute long time interval of the original 
data starting at 2:55:30 UT. The time resolution is 0.5 s and 
the frequency resolution is 1/256 Hz. This signature exhibits 
an extremely good signal to noise ratio. 
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The ideal detector has a threshold r with no false 
alarm rate {FA{t) — 0) and perfect detection efficiency 
{DE{t) = 1). This implies Nq = Nd{t). The optimal 
threshold r is the best compromise between detection ef- 
ficiency and false alarm rate. A threshold with zero false 
alarm rate and a small detection efficiency is not desir- 
able because it means that only rare events are detected 
(perfect chirp signals unaffected by noise). In this case 
the number of detected occurrences are affected by noise. 
The same is true for a perfect detection efficiency but a 
large false alarm rate where many detections result from 
signatures other than chirps. 

Therefore we investigate the properties of the detec- 
tion efficiency (DE) and the false alarm rate (FA) as a 
function of the threshold t for each detector. 



The standard deviation detector 



III. AUTOMATED DETECTORS AND THEIR 
PERFORMANCE 

The performance of an automated detector operating 
with a threshold t is measured by the Detection Effi- 
ciency DE{t) and by the False Alarm rate FA{t) . The 
number of sprites occurring Nq can be obtained from the 



The standard deviation detector detects a chirp when 
the signal intensity exceeds a given threshold. We calcu- 
late the signal intensity changes in the high pass filtered 
infrasound recording, by moving a window of length Ls 
through the data and evaluating at each time the stan- 
dard deviation of the data inside the window. We choose 
Ls = 256 data points. This corresponds to a 12.8 seconds 
long time interval: a time interval in the expected range 
of the average duration (T~12-15 s) of a sprite chirp sig- 
nature. The variation of the standard deviation is shown 
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in Fig. The horizontal line indicates the threshold 
T — 0.015 Pa., while the bottom diagonal line indicates 
the long term trend of the standard deviation. The in- 
tensity follows the day-night cycle of temperature, with 
a maximum around noon and a minimum around mid- 
night. In Fig. 13 the intensity slowly increases as the sun 
rise approaches at 4:00 UT. 

A detector based on the standard deviation should con- 
tinuously scale the threshold to take in account the ef- 
fect of change of intensity in time or particularly windy 
conditions. Here we use the standard deviation detec- 
tor for comparison with the spectral covariance detector 
and the CASSANDRA detector. Therefore we use a con- 
stant threshold for the entire duration (2 hours) of the 
infrasound recording. 

Fig. ^ shows the detection effieciency and the false 
alarm rate for different values of the threshold t. It 
is evident that the standard deviation detector has no 
good compromise between false alarm rate and detec- 
tion efficiency. For values of r < 0.017 we have a de- 
tection efficiency (DE > 0.6 or 60%) but a false alarm 
rate (FA > 0.75 or 75%). Raising the threshold lowers 
the false alarm rate to about 0.6 (60%) but the detection 
efficiency decreases to less than 0.3 (30%). These results 
for the standard deviation detector are not surprising: 
every chirp signature implies an increase in the signal to 
noise ratio but the inverse conjecture is not true. 
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FIG. 3: The variation of the standard deviation a of the high 
pass filtered signal in a 12.8 second long interval. The hor- 
izontal dashed line indicates a threshold of 0.015 Pa. The 
diagonal dotted line at the bottom shows the increase in in- 
tensity of the pressure fluctuations as sun rise approaches. 



B. The spectral covariance detector 

The spectral covariance detector rests on the appli- 
cation of the spectral covariance of the spectrogram to 
detect chirp signatures of increasing frequency (Fig. [2J|. 
In this section, we introduce the spectral covariance and 
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FIG. 4: Standard deviation detector: the Detection Efficiency 
(squares) and the False Alarm rate (circles) as a function of 
the threshold r. 

discuss the performance of the spectral covariance detec- 
tor for the infrasound chirp signatures of sprites. 



1. Spectral covariance 

The Spectrogram of a signal is its energy density 
E{t, /) in the time-frequency plane. A linear chirp of 
increasing frequency will appear as a diagonal line in the 
spectrogram (Fig. EJ. The inclination of the line with 
respect the horizontal axis is proportional to the rate v 
at which the frequency of the chirp changes in time. The 
spectral covariance uses the covariance between the en- 
ergy density E{t, /) relative to two different frequencies 
to detect the presence of a diagonal line in the spectro- 
gram. The covariance of time delay S relative to the 
frequencies / and / + fs is 

Cfj^ (5) =< E{t, f)E{t + SJ + fs) >t, (2) 

where < ... >t denotes the time average and /s>0 is the 
frequency shift. 

For a white noise signal, the computed spectrogram 
E{t^ f) is a randomly fluctuating function in both the ar- 
guments t and /. Thus, every delay d has the same prob- 
ability of maximizing the covariance Cfj^{5). The same 
holds true if a periodic component of fixed frequency fo 
is superimposed on the noise. In this case, if both / and 
f+fs are different from fo then E{t, f) and E{t-\-5, f+fs) 
are randomly fluctuating functions, if f=fo then E{t, fo) 
is constant but E{t + S, f + fs) is a randomly fluctuat- 
ing function and vice-versa. In the case of a linear chirp 
with initial frequency fi and final frequency // > fi the 
covariance Cfj^ (S) will have a maximum at Smax 7^ 
whenever the frequencies / and / -I- fs are in the range 
[/»://] or, equivalently, 

fs<ff~h- (3) 
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The value Smax is inversely proportional to the rate v at 
which the frequency of the chirp is changing: 



V 



(4) 



Finally, in the case of an impulsive signature (a straight 
vertical line) in the spectrogram like the one of the 
thunder produced by a lightning, we expect 5max=^ 
(?;—>• +00). 

To use the covariance of Eq. (0) for detecting linear 
chirps with a variable frequency range we need to elimi- 
nate its dependence on a particular frequency / without 
loosing the useful properties in chirp detection. Thus we 
average the covariance over all frequencies and define the 
spectral covariance of delay 5 and frequency shift fs as 



SCfS5) ^<CfjM >/, 



(5) 



where < ... >f denotes the frequency average. The delay 
Smax for which the spectral covariance of Eq. holds its 
maximum value, has exactly the same properties of the 
delay S„iax relative to the covariance of Eq. ((SJ. In Fig.[Sl 
we plot the covariance C/j^ and the spectral covariance 
SCf^ as a function of the delay 6 for the 2 minute long 
spectrogram of Fig.|21 The chirp of Fig.|21produces large 
values of the covariance for frequencies / between 3 and 4 
Hz and delays S between 5 and 10 seconds. The spectral 
covariance has its maximum for 6max—6 seconds. This 
value of 6max corresponds (Eq. ^ to a rate of frequency 
change w~0.13 Hz/s and an inclination of the chirp sig- 
nature in the spectrogram of ~58° (Fig.[5J). 




FIG. 5: (a) The covariance C/,/^ (S) as a function of the delay 
S (in seconds) for the spectrogram of Fig. 2. Here /s = 10/256 
Hz (~0.8 Hz), (b) The spectral covariance SCf^ as a function 
of the delay 5 (in seconds). 



the location t is examined and the delay Smax(t) evalu- 
ated. Then the interval is shifted and centered to the new 
location t + u and the correspondent Smax {t + u) is eval- 
uated. Consecutive intervals containing a linear chirp of 
increasing frequency will result in consecutive equal val- 
ues of the delay 6max{t) ^ 0. For a linear chirp of time 
duration T, the number ncEV of consecutive equal values 
of the delay Smax is in the range 



T ^ I 

— < ncEV ^ -■ 
u u 



(6) 



In real cases a chirp signature in the spectrogram will 
produce similar consecutive values of the delay 6max{t)- 
The automated detection algorithm measures the disper- 
sion of ncEV consecutive values of Smax{t) around their 
mean. A detection will be reported if the relative disper- 
sion (the ratio between the standard deviation and the 
average) is below a given threshold r. 

There are some numerical limitations in the evalua- 
tion of the spectral covariance SC which must be con- 
sidered. (1) Numerically, one evaluates the energy con- 
tained in a box of dimension AtxAf centered in the 
location {t,f) of the time-frequency plane 0|. The nu- 
merical energy density E^uMit, f) is obtained dividing 
the energy by the dimension of the box. The locations 
{t, /) for which the numerical energy density ENUM{t, f) 
is computed belong to a grid of steps St in the time do- 
main and Sf in the frequency domain. The intervals 
St and Sf are the time and frequency resolution of the 
spectrogram, while the intervals At and A/ are the time 
and frequency localization of the spectrogram. Typically, 
At>6t and Af>5f. These numerical limitations impose 
a lower bound on the frequency shift fg 



fs > A/. 



(7) 



If the condition of Eq. ^ is not satisfied, the numerical 
energy densities ENUM{t,f) and EMUM{t,f + fs) used 
in the evaluation of the spectral covariance (Eq. (S)) refer 
to two overlapping intervals of frequencies. In this case 
Eq. may not be satisfied. (2) The duration / imposes 
a limitation on the values of delays S for which a sta- 
tistically meaningful numerical evaluation of the spectral 
covariance, is possible. If Sn / is the maximum delay 
used in evaluating the spectral covariance, only linear 
chirps with a rate of frequency change of 



(8) 



can be detected. 



3. Application to infrasound from sprites 



2. Ctiirp detection 

The spectral covariance detector operates as follows. 
An interval of the spectrogram of duration / centered in 



A visual inspection of the chirp signatures of the 9 
sprites present in the infrasound recording examined, 
shows that the difference between the final and the ini- 
tial frequency of the chirps is «2-2.5 Hz. We evaluate the 
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spectrogram with a frequency resolution of 1/256 Hz and 
we set the frequency shift fs to 10/256 Hz (0.78125 Hz). 
This value of the frequency shift satisfies the inequali- 
ties of Eqs. lEJl and Q. We select 7=30 s and 6r=20 s. 
With this value of Sr only chirps with a rate of frequency 
change v greater than 0.039 Hz/s can be detected. These 
chirps will produce signatures of inclination greater than 
27° in a 2 minutes long display of the spectrogram. The 
shift u of two consecutive intervals of duration / is chosen 
to be equal to 3.2 seconds and consequently (Eq. we 
set ricEV='^- Finally we want to exclude the possibility 
to detect impulsive signature from lightning. This signa- 
tures in theory should produce a Smax=0 s, corresponding 
to vertical lines (inclination of 90°) in the spectrogram. 
In practice, however, it is better to exclude 6max too close 
to zero. We consider only those delays Smax>i-8 s cor- 
responding to an inclination of ~80° in a 2 minute long 
display of the spectrogram. 

In Fig. 1^1 we plot the sequence of delay S„iax{t) for 
the same 2 minute long time interval of the infrasound 
recording used for Fig. |21 The presence of the chirp co- 
incide with consecutive almost equal values of the delay 
Smax{t) (around t~60 s). Before and after the chirp sig- 
nature a small number of consecutive similar values of 
Smax{t) (around t—20 s and i=100 s) and some isolated 
fluctuating values. In Fig. |7| we plot the detection effi- 
ciency and the false alarm rate for different values of the 
threshold r. As for the standard deviation detector there 
is no good compromise between the false alarm rate and 
the detection efficiency. For values of r > 0.1 the spectral 
covariance detector has a detection efficiency of almost 
80%, but a false alarm rate superior to 90%. Lowering 
the threshold results in a false alarm rate slightly below 
80%, but in a drop of the detection efficiency from ~80% 
to ~20%. This behavior of the detection efficiency and of 
the false alarm rate has two causes. (1) "Spurious" sig- 
natures (not from sprite) produce a sequence of values of 
Smax with a small (<0.05 or 5%) relative dispersion. (2) 
Sprites signatures may produce sequence of delay S„iax 
with an large (>0.1 or 10%) relative dispersion. 
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FIG. 6: The delay Smax (in seconds) for a interval of duration 
I of the spectrogram centered at the location t (in seconds). 
This plot refers to the same 2 minute long interval of the 
original data used for Fig. 2. The dashed line at the bottom 
indicates the delays Smax=i-8 s expected for a signature with 
a inclination of ~80° in a 2 minutes long spectrogram. 
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FIG. 7: Spectral covariance detector: the Detection Efficiency 
(squares) and the False Alarm rate (circles) as a function of 
the threshold r. 



C. The CASSANDRA detector 



The detector described in this section derives its name 
from the CASSANDRA analysis [ulU^ : an apphcation 
of diffusion entropy analysis Il3l Il4j to non stationary 
time series In the following we briefly discuss the 

details of the diffusion entropy analysis and the changes 
to the original [H [l^ formulation of the CASSANDRA 
analysis necessary for the detection of chirps. We then 
show the results of the application of the CASSANDRA 
detector to the infrasound chirp signature of sprites. 



1. Diffusion entropy analysis 

The application of the diffusion entropy analysis to a 
time series is made of two steps. Step 1: use the time 
series to create a diffusion process. Step 2: monitor, 
as diffusion takes place, the entropy of the probability 
density function (pdf) describing the diffusion process. 
The behavior of the entropy is indicative of the statistical 
properties of the time series analyzed. 

The first step of the diffusion entropy analysis is com- 
puting all the possible sums of any n consecutive terms 
of the time series {^j} of length L, namely: 
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k+n-l 

Xk,n = E ^3 (9) 
j=k 

k = I, L ^ n + I. and n~l,...,L. 

This procedure describes a diffusion process if we con- 
sider the sequence S,j as the sequence of fluctuations of 
a diffusing trajectory and n as the time for which the 
diffusion process has taken place. Consequently, each of 
the L—n+1 values Xk,n can be thought as the position 
of a diffusion trajectory after a time n starting from the 
location at n—0. The second step is computing the pdf 
p{x, n) of finding a trajectory in the location x after a 
time n and its diffusion entropy 

S{n) = — J p{x,n)\n p{x,n)dx. (10) 

The numerical evaluation of S{n) is done by dividing at 
each time n the diffusion space in cells of equal size A(n) 
centered around the location Xj. The size A(n) must be 
small enough for the pdf p{x, n) to be constant inside the 
cell. In this case, 

S{n) ~ — ^^^(xj , n) hip{xj,n) + In A(n) (11) 
j 

with 

p{xj,n) = p{xj,n)A{n). 

The choice of having a temporal dependence on the cell 
size A is due to the necessity of satisfying the condition 
p{xj, n) constant inside a cell and the probability p(a;j, n) 
being large enough for a meaningful statistical evaluation 
with L—n+1 trajectories at time n (Eq. (SJ. A choice of 
a small fixed size cell will satisfy the former condition 
but not the latter when n increases and the diffusion 
trajectories explore larger intervals of the diffusion space. 

A time series {^j} of random uncorrelated numbers 
drawn from a distribution of finite variance generates a 
diffusion process that rapidly becomes Brownian . As 
a consequence, the diffusion entropy S{n) will rapidly 
approach a regime of linear increase on a logarithmic time 
scale with slope 0.5 j^O]- 

The addition to the random time series of a pe- 
riodic component {^*^} with a periodicity of M data 
points, has the effect of "bending" the diffusion entropy: 
the periodic component increases the value of S(n) at 
times n that arc not multiples of M and has no effect 
on the value of S{n) at times n multiples of M. As a 
consequence, the times n multiple of the period M are 
now points of local minima for the diffusion entropy S{n). 
The "bending" effect is shown clearly in Fig. |S1 When 
a periodic component is added to a random noise 



fluctuation {^j} we can write the sum Xk.n of Eq. ^ as 

k+n-l 

j=k 

k-\-n—l k-\-n—l 

= E 0+ E ef = :^rr+<r*^ (12) 

j=k j=k 

k = 1, L — n + 1. and n~l,...,L. 

The addition of the periodic component {i^*^} has the 
effect of shifting the position x^'J^*^ at time n of the A:-th 
diffusion trajectory by the quantity If n is not 

a multiple of the period M, the shift ^p^™'^^'^ depends on 
the index k. Therefore different diffusion trajectories are 
shifted by a different quantity. This results in a bigger 
"spreading" of the diffusion trajectories and therefore in 
a larger value of the diffusion entropy S{n). When n is 
a multiple of the period Af , the shift x^^"° ^'^ is inde- 
pendent from the index k. In this case all the diffusion 
trajectories are shifted by the same quantity and the dif- 
fusion entropy S{n) does not change. 

Moreover, Fig. |S1 shows how the "bending" of the dif- 
fusion entropy S{n) becomes smaller and smaller as the 
time n increases. The standard deviation of the sum 
^k°n"^ (Eq. I12|) of n consecutive terms of the random 
fluctuations {^j} increases as n increases. The standard 
deviation of the sums ^^^"'^"'^^'^ (Eq. I12|) of n consecutive 

terms of the periodic component {C/^} is limited and it 
is periodic of period M: it vanishes whenever n is a mul- 
tiple of the period M, it increases, reaches a maximum 
and then decreases in between two consecutive periods. 
Therefore the contribution to the diffusion entropy S{n) 
of the periodic component becomes smaller and smaller 
compared to that of the random values {^j } and so does 
the "bending" effect. 

2. CASSANDRA analysis 

The CASSANDRA analysis is the application of the 
diffusion entropy analysis to smaller intervals of a time 
series such that the differences between the results in 
each interval reflect the statistical changes occurring in 
the time series itself. To compare the differences between 
the diffusion entropy of different intervals of a time series, 
the authors of ^11, 2^ use the following quantity 

= ]^ E [^*(") - ("^^ (1) + Hn))] ■ (13) 

n=l 

In Ea. (|13|l . t is the position in the time series where the 
small interval of length Lg is centered, 5*4 is the diffusion 
entropy relative to this interval of the time series and N 
is the maximum time for which the evaluation of St (n) is 
statistically meaningful. Cl^ (t) is the difference between 
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FIG. 8: The difFusion entropy S{n) as function of time n 
of gaussian noise (line and white squares) and of gaussian 
noise plus a periodic component (line and black squares) . The 
periodic component has a periodicity of Af of 8 data points. 
The full line represents the increase of S{n) with a slope of 
0.5 on a logarithmic time scale. 



the "local" diffusion entropy S't(n) and one that, starting 
from the same value at n = 1, increase with a slope of 0.5 
on a logarithmic time scale. Therefore, Cl, {t) is an indi- 
cation of how different the diffusion process generated by 
the data in the small interval centered at t is from Brown- 
ian diffusion. The quantity Cl^ (t) of Eq. H13|l is useful in 
detecting increases of the diffusion entropy St{n) with an 
average slope smaller or larger than 0.5 as a result of the 
local correlation properties. But in our case, we want to 
detect chirps in background noise. Intervals containing 
a chirp will result in a "bended" diffusion entropy with 
the times n of the local minima of St{n) depending on 
the chirps instantaneous frequency. For this reason we 
evaluate, instead of the quantity defined in Eq. ((T^ . the 
cumulative slope change CSClS'^) and the time fiL^{t) 
of the first local minima of St{n). 

The cumulative slope change for a interval of length 
Ls centered at the position t of the time series is defined 
as 



CSCLAt) 



1 



N-2 

E 



t{n+l)-et{7i)\ 



TV - 2 ^ ln(n + 2) - \n{n) ' 



where 



tin) 



St{n + l)- Stjn) 
\n(n + 1) — ln(n) 



(14) 



(15) 



St{n) and N of Eq. (|14|l are the same quantities as in 
Eq. ((T^ . while 9t{n) of Eq. is the slope of the hnc 
connecting two consecutive values of the diffusion entropy 
St{n) when plotted on a logarithmic time scale. The cu- 
mulative slope change of Eq. H14|) is the weighted sum of 
the absolute value of the difference between two consec- 
utive slopes. The weights are the intervals of time (on 



a logarithmic time scale) during which the slope differ- 
ence is evaluated. Slope changes happening at later times 
are weighted more. This compensates the fact that the 
"bending" of the diffusion entropy becomes smaller as n 
increases. The cumulative slope change is able to dis- 
tinguish intervals of the time series where the diffusion 
entropy St{n) is "bended" from those where it is not. 
The time of the first local minima (t) is the first time 
n for which the conditions 9t{n — l) < and 9t{n) > are 
satisfied. If these condition are not met for any n < N 
then we set nL^t) = N. Therefore for an interval of 
length Ls centered at the position t which contains only 
noise rii^ (t) = N jsj . For a sequence of intervals with a 
periodic component of fixed period M < N, fiL^t) = M. 
Finally for a sequence of intervals containing a chirp sig- 
nature fiL^t) will change accordingly to the local chirp 
frequency. 



3. Chirp detection 

For the purpose of chirps detection we need to use 
the cumulative slope change CSCls (t) together with the 
time of first local minima n^^ (t) . . We define the chirp 
cumulative slope change cCSCL,{t) as 



Oif • 



cCSCLsit) = 



flLsit)-- 
flLsit)-- 



--N 



t e J 



(16) 



^CSCL^t) otherwise 



The set J is the set of the all the locations t for which 
consecutive constant values oi nLs{t) = M < N are en- 
countered. The chirp cumulative slope change of Eq. 1)16(1 
vanishes for intervals containing only noise or a periodic 
component of fixed period M. The CASSANDRA de- 
tector detects a chirp when the chirp cumulative slope 
change cCSCl^ (t) exceeds a given threshold r. 



4- Application to infrasound from sprite 

As for the standard deviation detector we choose Lg = 
256 corresponding to a 12.8 seconds long time interval. 
Moreover, we dichotomize the signal such that every data 
point above the average is -1-1 and every data point be- 
low the average is —1. This drastic procedure has three 
advantages. (1) The night-day cycle of the intensity in 
the infrasound recording (Fig. |2Jl is eliminated by the 
dichotomization. Thus, a unique threshold independent 
from the time of day can be set for the purpose of auto- 
mated detection. (2) The problem of choosing an appro- 
priate value of A(n) for different times n in Eq. I|ll|) is 
simplified: a unitary cell will be used for all the values of 
n explored. (3) The signal to noise ratio is "preserved" 
by the dichotomization: a very intense sprite signature 
will result in an almost perfect periodic pattern with de- 
creasing frequency, while a weak signature will be difficult 
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to recognize because the noise will randomly affect the 
pattern. 

In Fig. El we plot the chirp cumulative slope change 
cCSCL^{t), the cumulative slope change CSC{t) (for 
clarity the cumulative slope change has been moved down 
with respect the chirp cumulative slope change) and the 
time of the first local minima n(t) for the same insert of 
2 minutes used for Fig. [21 We clearly see that the chirp 
signature produces a big value of the cumulative slope 
change and that the time of first local minima n^^ detects 
its change in frequency passing with continuity from 7 to 
4. The chirp cumulative slope change of Eq. Hlt)|) reduces 
the possibility of a false alarm annulling the cumulative 
slope change in the case of noise or periodic component 
with fixed periodicity., 

Finally, in Fig. EH we plot the the detection efficiency 
and the false alarm rate for different values of the thresh- 
old T. The CASSANDRA detector has a good compro- 
mise between the false alarm rate and the detection effi- 
ciency. For a threshold value of 0.8 the false alarm rate 
is null and the detection efficiency is about 66%. 





0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 



FIG. 10: CASSANDRA detector: the Detection Efficiency 
(squares) and the False Alarm rate (circles) as a function of 
the threshold r. 



would expect to get a better false alarm rate lowering the 
threshold, but this is not the case (panel (b)). Finally we 
see how raising the threshold in the CASSANDRA detec- 
tor improves the false alarm rate without compromising 
the detection efficiency (panel (c)). 



FIG. 9: The cumulative slope change CSCls (squares), the 
time h (full line) of the first local minima of the diffusion 
entropy S*(n) and the chirp cumulative slope change cCSCls 
(dotted line) as a function of the position t (in seconds) . The 
plot refer to the same 2 minute long interval of the original 
data used for Fig. 2. The cumulative slope change has been 
translated down for clarity. The horizontal line at the bottom 
indicates a threshold of 0.8. 



IV. CONCLUSION 



The plots of the detection efficiency and false alarm 
rate as a function of the threshold r (Figs. 01 Hand EH) 
indicate that the CASSANDRA detector is the one with 
the best trade off between detection efficiency and false 
alarm rate. This is confirmed by the plot of Fig. ^2 We 
see how raising the threshold r in the case of the stan- 
dard deviation detector (panel (a)) does not lower the 
false alarm rate. For the spectral covariance detector we 
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FIG. 11: Comparison between the three different detectors 
used, (a) The standard deviation detector, (b) The spectral 
covariance detector, (c) The CASSANDRA detector. In all 
cases the the dashed vertical line indicate the time of the 
sprite signatures, while the horizontal dotted lines indicates 
a threshold. 
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